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ABSTRACT 

In the X-ray spectra of most X-ray dim isolated neutron stars (XDINSs) ab- 
sorption features with equivalent widths (EWs) of 50 - 200 eV are observed. 
These features are usually connected with the proton cyclotron line, but their 
nature is not well known yet. We theoretically investigate different models to 
explain these absorption features and compare their properties with the observa- 
tions for a better understanding of the radiation properties of magnetizedneutron 
star surfaces. Based on these models, we create a fast and flexible code to fit 
observed spectra of isolated neutron stars. We consider various theoretical mod- 
els for the magnetized neutron star surface: naked condensed iron surfaces and 
partially ionized hydrogen model atmospheres, including semi-infinite and thin 
atmospheres above a condensed surface. Spectra of condensed iron surfaces are 
represented by a simple analytical approximation. The condensed surface radia- 
tion properties are considered as the inner atmosphere boundary condition for the 
thin atmosphere. The properties of the absorption features (especially equivalent 
widths) and the angular distributions of the emergent radiation are described for 
all models. A code for computing light curves and integral emergent spectra of 
magnetized neutron stars is developed. We assume a dipole surface magnetic 
field distribution with a possible toroidal component and corresponding temper- 
ature distribution. A model with two uniform hot spots at the magnetic poles 
can also be employed. Light curves and spectra of highly magnetized neutron 
stars with parameters typical for XDINSs are computed using different surface 
temperature distributions and various local surface models. Spectra of magne- 
tized model atmospheres are approximated by diluted blackbody spectra with 
one or two Gaussian lines having parameters, which allow us to describe the 
model absorption features. The EWs of the absorption features in the integral 
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spectra cannot exceed significantly 100 eV for a magnetic atmosphere or a naked 
condensed surface as local surface models. A thin atmosphere above a condensed 
surface can have an absorption feature whose EW exceeds 200 eV in the inte- 
grated spectrum. If the toroidal component of the magnetic field on the neutron 
star atmosphere is 3 - 7 times larger than the poloidal component, the absorp- 
tion feature in the integral spectrum is too wide and shallow to be detectable. To 
explain the prominent absorption features in the soft X-ray spectra of XDINSs a 
thin atmosphere above the condensed surface can be invoked, whereas a strong 
toroidal magnetic field component on the XDINS surfaces can be excluded. 

Subject headings: radiative transfer - scattering - methods: numerical - (stars:) 
neutron stars - stars: atmospheres - X-rays: stars 
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Introduction 



X-ray dim isolated neutron star s (XDINSs) represe nt a new class of X-ray sources 



established by ROSAT observations ( iWalter et al. 



19961). At present, se ven such objects 



Haberl 



20071 ). All of these 



have been found ("The Magnificent Seven", see review by 
objects have ages ~ 10^ yr and thermal-like spectra with effective temperatures Tes ~ 10^ 
K. The radiation of XDIN Ss is pulsed with a mplitudes from 1% ( RX J1856.5— 3754, 



Tiengo fc Mereghetti 



20071 ) to 18% (RBS 1223, ,Schwope et al 



20051) 



Abso rption 



XDIN Ss (iHaberl 



eatur es in the 0.2-0.8 keV range are detected in the spe ctra of almost a ll 



20071 ). In t hree objects the presence of two (RBS 122 3. ^Schwope et al. 



Haberl 



2003) or three (RX J1605.3+3249, 



Haberl 



20071 ) 



20071 . and RX J0806.4-4123, 
absorption lines is claimed . Equ ivalent widths of these features are about 30-100 eV 



(IHaberl et al. 



eV 



200' 



Schwope et al 



Haberl 



20071 ) with the most prominent absorption in RBS 1223 (~ 200 
20071 ). If these features are interpreted as proton cyclotron lines, they 
correspond to magnetic fields B ~ 0.4 - 1.6 xlO^^ G on the XDINS s surfaces, but if they are 



interp reted as electron cyclotron lines, then i? ~ 2 - 8 xlO G (see 



20071 ). 



van Kerkwijk fc Kaplan 



The optical counterparts of s ome XDINSs have been fo und (see reviews 



Mignani et al. 



2007 



Mignani 



2009 



and references therein, and 



by 



Eisenbeiss et al. 



20101 1. 



It is very intriguing that their optical/ultravio 



blackbody extrapo 



2003 



Motch et al 



ation of the X-ray spectra (iBurwitz et al. 



et fluxes are a few t imes l arger than the 



2001 



2003 



Kaplan et al. 



20031). Two explanations for this excess have been suggested. First, it 



can be caused by non-uniform temperature distributions across XDINS surfaces. Second, 
the emergent spectral energy distribution of the surfaces can differ from the blackbody. For 
example, using non-magnetic model atmosphere sp e ctra gives an ey en larger optical excess 



compared to the observed one (jPavlov et al. 



19961) 



Ho et al. 



(120071 ) explained the optical 
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excess in RX J1856.5— 3754 by a single magnetic model atmosphere of finite thickness above 
a condensed surface. Both of these mechanisms can be important for XDINSs. 

The XDINSs are nearb y objects, and distances to some of them have been estimated 



by parallax measurements (iKaplan et al. 



helps to evaluate their radii ( iTriimper et al. 



2002 



van Kerkwijk Sz Kaplan 



20071 ). This 



20041 ). yielding useful information on the 



equation of state (EOS) for the neutron-star core, which is i mportant for the studies o f 



particle interaction s at supranuclear densities and pressures flLattimer &: Prakash 



2007 



Haensel et al. 



20071 ) as well as for computatio ns of templates of g ravitational wave signals 



which arise during neutron stars merging (e.g. 



Baiotti et al. 



20081) 



The theoretical emergent spectra of the XDINS surfaces are also necessary for the 
evaluation of neutron star radii. The XDINS surface laye rs can ei ther be condensed or 



haye plas r na atmospheres ( 



Lai 



2001 



Medin fc Lai 



Ruderman 



1971 



Romani 



1987. see also 



Lai fc Salpeter 



19971 : 



20071 ). depending on surface temperature and magnetic field 



strength. The radiation spectra of magnetized cond ensed surfaces we r e calculated by ma ny 



authors in the frarnework of different assumptions 



Perez- Azorm et al. 



2005 



van Adelsberg et al. 



(IBrinkmann 



1980 



TuroUa et al. 



2004 



20051 ). Investigations of the neu tron-star 



atmospheres are even more extended, from the first relatively simple models (iRomani 



1987 



2002 



Shibanov 



Ho fc Lai 



et al 



19921) to m ore sophisticated models, including comple te (jLai fc Ho 



20031 ) or partial flLai fc Ho 



2003 



yan Adelsberg &: Lai 



conversion due to the vac uum polarization effect (Pavloy fc Gnedin 



ionized hydrogen models ( jPotekhin et al. 



2004; 



Ho et al 



20061) photon mode 



19841 ). and partially 



20081 ). Mid-Z element partially 



ionized atmospheres for magnetized neutron stars (however without detailed account 
of atomic motion effects) have also been modeled and used for the absorption features 



interpretation (see iMori &: Ho 



20071 and references therein) . 



Absorption features in the soft X-ray spectra of XDINSs may be the Rosetta stone for 
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understanding of the physical nature of their surfaces. For this purpose it is necessary to fit 
the observed spectra by theoretical spectra of model atmospheres and condensed surfaces. 
A correct model should allow one to describe the observed absorption features, as well as 
the overall spectral energy distribution, including the optical/UV fluxes, and the pulse 
profiles. For such a comparison an integration of local spectra over the neutron star surface 
should be performed with relativistic effects taken into account. Therefore, the distribution 
of the magnetic field and the effective temperature over the neutron star surface has to 



be found also. Computations of integral 
were carried out by many authors (e.g.. 



spectra from loca . 



Zavlin et al 



1995 



atmosphere eraergent spectra 



Zane fc TuroUa 



2006 



Ho et al. 



20081 ). Some a uthors have also usee 



X-ray spectra (jZane fc TuroUa 
therein). 



2006 



model at mospheres to fit th e observed is olated XDINS 



Ho 



20071 . see also review by 



Zavlin 



2009 



and references 



Usually a simple dipole magnetic field distri bution over the neutron star sur: 



corresponding relativistic correction is assumed (jPavlov &: Zavlin 



200G 



Ho et al. 



ace with 



20081) 



If this surface dipole field corresponds to a global dipole field of the star, one can find 
a surface effective temperature distribut ion depending on the magnetic field strength 



and 



: he inner neutron-star temperature (IGreenstein fc Hartke 



2001 



1983 



Potekhin fc Yakovlev 



Potekhin et al. 



20031 ). The surface temperature distribution is defined by the local 
inclination angle of the magnetic field to the surface normal, which determines the value 
of the local radial thermal conductivity in the neutron star crust. Bright pole spots are 
too large to explain the observed XDINSs' pulsed fractions if the dipole field model a nd 



the blackbody approximation are used for the emergent spectra (IGeppert et al. 



20061). 



Quadrupole or more comp 
have been suggested (e.g.. 



icated magnetic field configurations in neutron star crusts 



Ruderman 



1991 



Zane fc TuroUa 



20061 ). which may provide 



an explanation for the high pulsed fraction, although an e qually likely explana tion 



possibly, may be provided by models of polar caps heating. 



Geppert et al. 



(120041 ) assumed 
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that the magnetic 



Perez- Azorm et al 



field is concentrated in the crust only. iGeppert et al.l ( 120061 ) and 



(l2006al ) propose that the magnetic field in the crust additionally has a 



strong toroidal component. In both cases t he polar hot spots are sma ll and can explain the 



observed XDINS pulsations. In particular, 



Perez- Azorm et al. 



( l2006bl ) used the model with 



the toroidal field to interpret all the observed features of RX J0720— 3125 on the base of a 
condensed surface emission. 

At present time, most of the observed XDINSs spectra were fitted by a simple blackbody 
model with one or more Gaussian lines. The detailed modeling of observed spectra by using 
accurate model atmosphere or condensed surface spectra is computationally very expensive, 
particularly if the temperature and magnetic field distributions over the stellar surface 
are considered. For example, the spectrum of RX J1856.5— 3754 was described by a single 



model atmosphere ( iHo et al. 



2QQ71) . and the 



by using only four latitude points (IHo 



pulsed fraction of the same object was modeled 



20071 ). This is one possible way for fitting spectra of 
XDINSs: Create an extensive set of theoretical XDINSs spectra with (inevitably) limited 
variation of stellar pa rameters as we ll as temperature and magnetic field distributions over 



the surface (see, e.g. 



Ho et al. 



20081 ). We suggest an alternative way. The local spectra 
of the neutron star together with temperature and magnetic field distributions can be 
fitted by simple analytical functions. Then an integral spectrum of the neutron star with 
a sufficiently detailed latitude set can be computed rather quickly. This approach is more 
flexible and faster, but the results will be more approximate. However, in this way one can 
first constrain the XDINS parameters, and then the more accurate, expensive calculations 
can be performed to refine the analysis. Here we develop this approximate approach and 
its application for fitting XDINSs spectra will be presented in a different paper. Some 
preliminary conclusions about the possible nature of the XDINSs emitting surface will be 



made on the base of test calculations. 
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In this work we analyze the properties of emergent spectra of magnetized model 
atmospheres, condensed surfaces, and thin magnetized model atmospheres above condensed 
surfaces, paying special attention to the equivalent widths of the absorption features and 
the angular distribution of the emergent radiation. Spectra of naked condensed surfaces 
are represented using a simple analytical approximation. We also present test calculations 
of the integral spectra and light curves of the magnetized neutron star models, using a 
simple representation of the considered local model spectra by analytical functions. These 
calculations allow us to formulate the qualitative characteristics of the local spectra that 
are necessary to fit the observed pulsed fractions of XDINS light curves and the absorption 
features of XDINS spectra. These simple computations allow us also to choose the model of 
the radiative neutron star surface, which will be used later for more extensive and accurate 
computations of the integral spectra and the light curves of XDINSs. 



2. Local models of a neutron star surface 



A magnetized neutron star surface can be condensed or can have a plasma atmosphere. 
It depends mainly on chemical composition, tem perature, and magne t ic field strength, but 



necessary conditions are not established exactlv (ILai fc 



Lai 



2001 



Potekhin &: Chabrier 



2004 



Medin fc Lai 



Salpeter 



1997 



Potekhin et al. 



1999 



20071 ). In this section we shortly describe 



theoretical radiation properties of the condensed surfaces and the model atmospheres of 
magnetized XDINSs. The main attention is concentrated on the equivalent widths of the 
absorption features and the angular distributions of the emergent spectra, which are also 
very important for the integral spectra computations. The model atmospheres are discussed 
in more detail. 
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2.1. Radiation from the condensed surface of magnetized neutron stars 



The description of magne t ized c ondensed surface radiation propert ies is based on the 



resu 



2005 



ts of 



van Adelsberg et al. 



(|2005[ ) (see their Figs. 2-6; the results of 



Perez- Azorin et al 



are in good agreement). They were obtaine d using a free- i on ap proximation, which 



can be not exactly accurate. An earlier work by iTurolla et al.l ( 120041 ) was based on the 
approximation of fixed (non-moving) ions. The real radiation properties of a condensed 
magnetized surface can b e intermediate between these limits (see the discussion in 



van Adelsberg et al 



20051 ) 



In the free-ion approximation, there exist two broad absorption features in the 
spectrum. The first one lies between the ion cyclotron energy 

Z 



Ec,i = h ZeB/m^c ^ 0.0635 B13 ( - ) keV, 



(1) 



where A is the mass number of the ion (f» 1 for H and ~ 56 for Fe), Z is the ion charge, mj 
is the ion mass, -B13 = G, and some boundary energy Eq 



Ec^ E,-, + ElJE, 



(2) 



where i^p.e is the electron plasma energy 



Ep^e = ^ 



rrip 



1/2 



0.0288 



A 



1/2 



P 



1 g cm 



1/2 



keV, 



and Ec^e is the electron cyclotron energy 



(3) 



Ec,e = heB/m^c ^ 115.8 B13 keV. 



(4) 



Here ne is the electron number density, p is the mass density, and me is the electron mass. 

Only one radiation mode can propagate between i?c,i and Eq, therefore the radiation 
emissivity is reduced to approximately 1/2 in this range. The condensed surface density 
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depends on magnetic field and chemical composition (e.g.. lLaill200ll ) 

Ps ^ 8.9 X 10^ T]' A Z -^'^ B^l^^ g cm-^ 



(5) 



where ri' is a coefficient ranging probably between 0.5 and 1 (see discussion in 

herefore, the plasma energy at the surface can be expressed 



van Adelsberff et al. 


200 


5). T 


as ( 


van Adelsberg et al. 


2005) 



2.7„'l/2 ^1/5 R 3/5 



(6) 



We conclude that it is possible to estimate the equivalent width of this absorption feature 
as: 

EW ^ 0.5-^ 117 r/' (Z/26)2/^ 5^^^ eV. (7) 
In all calculations below we used t]' = 1. 



The second feature is at the plasma energy Ep^, which is not interesting for us. Indeed, 
the observed absorption lines lie at < 1 keV, which would correspond to the plasma 
energy at i? < 6.4 x 10^^ {Z/26)~^^^ G. Meanwhile, the critical temperature for the plasma 
ph ase transition leading to the for mation of the condensed surface is roug hly 7.5 X IO^eJ/^ 



K (IPotekhin et al. 



1999 



Lai 



20011 ). Noting that the observed temperatures of the XDINSs 
are about 10^ K, we conclude that the absorption feature at the plasma energy in the 
condensed surface radiation cannot be less than 1 keV in these stars. Therefore, the 
observed absorption features in XDINSs cannot be associated with the absorption feature 
at the plasma energy. However, the first absorption feature (between Ec;i and Eq) can be 
prominent in their spectra. 



The angular distribution o 



results of 



van Adelsberg et al. 



the condensed surface emissivity can be estimated from the 



(120051) 



R) 5p 
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where Be is Planck function, 



1-R = (1 - 0.27(1 - cos af) (1 - 0.36(1 - cos ^Y) 
for E > Ec^ , 

l-i? = 0.3 + 0.2cos$ for E^;, < E < Ec„ 
1-R = (0.5 + 0.25 cos $) for E < Ec,i, 



and 



E, 



Ci 



Ec (1 + 3(1 - cos a)^/^) 



(9) 

(10) 
(11) 

(12) 



Here a is the angle between radiation propagation direction and the surface normal, and $ 
is the angle between the surface normal and the magnetic field lines. This approximation 
is presented graphica lly in Fig. [1] together with accurate calculations, performed by the 



method described by |yan Adelsberg et al. 



(120051 ) in the free ions approximation. 



Here we neglect the angular dependence on the azimuthal angle between radiation 
propagation plane and the normal- magnetic field line plane for the inclined magnetic 



field (see Fig. 6 in 



yan Adelsberg et al. 



20051 ) and the absorption feature at E^^- We also 



ignore the damping effects, which smooth the transitions between different bands. That 
approximation leads to sharp boundaries at Ec^i and Eq^ in the computed spectra (Figs. 
[3l m [6] - [HI [13]), which can be smoother in more detailed computations (see, however, next 
paragraph). In our approximation we ignore the azimuthal dependence of R in the case of 
an inclined magnetic field. At some azimuthal angles the considered absorption feature can 



be less significant (see Fig. 4 in 



van Adelsberg et al. 



20051 ). Therefore, our calculations with 



this approximation will give an upper limit for the strength of the absorption feature. This 
effect must be less significant in the case of a thin model atmosphere above the condensed 
surface due to significant atmospheric optical depth in the X-mode around the proton 
cyclotron resonance (see Fig. [5]). 
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Our approximation appears rough, but the more accurate calculations of R have 
significant uncertainties, too. First of all, these calculations were performed in two different 
approximations: a complete neglect of ion motions, and free ions. The real picture is 
in some sense intermediate. Ions in the crystal grid (or in liquid) are not free, but can 
move. These two limits g ive very different spectra at energies below i?c,i (see Figs. 2-6 in 



van Adelsberg et al. 



20051 ). The approximately known factor rj', which determines the value 



of Eq and the width of the absorption feature, was already mentioned above. Moreover, 
the accurate calculations are only applicable for a perfectly smooth condensed surface. 
If the surface has some rou ghness, the emitted spectrum has to be close to a blackbody 



(Ivan Adelsberg et al. 



20051 ) . Our approximation corresponds to the free-ions limit and gives 
radiation spectra, that reflect the main features of condensed surface spectra, especially if 
we take into account all uncertainties mentioned above. 



2.2. Model atmospheres of magnetized neutron stars 



Model atmospheres for magnetized neutron stars were computed by many scientiflc 
groups, but we use our own c alculations for the descri ption of their radiation spectra. We 



Suleimanov et al. 



(120091 ). where details of the computational 



employ the code presented by 
method can be found. 

Here we assume that the magnetic fleld strengths of the XDINSs are about a few xlO^^ 
G and consider a pure hydrogen atmosphere. The eff ective temperatures are ~ 10^ K. At 



1999|) andt 



ne vacuum 



these values of B and T hydrogen is partially ionized (iPotekhin et al. 
polar i zation effect together with the partial mode conversion can be signiflcant (iHo fc Lai 



2003 



van Adelsberg fc Lai 



20061 ). Therefore the model atmospheres presented below 



include these effects. The magnetic fleld in presented models is assumed to be homogeneous 
and normal to the surface. In general, it is necessary to study model atmospheres with 
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Fig. 1. — Approximation of dimensionless emissivity as a function of photon energy E for 
the case of a condensed iron surface at i? = 10^'^ G. Top panel: The magnetic field is normal 
to surface. The different curves correspond to different angles a between emergent photon 
direction and surface normal. Bottom panel: The emergent photon direction is inclined 45° 
to the surface. The different curves correspond to different angles $ between magnetic field 
lines and surface normal. 
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inclin ed magnetic field for the modeling of integral spectra of neutron stars (iHo et al 



20081 ). But here we consider general features of the model spectra, and therefore we neglect 
the distinctions between spectra of the models with different field inclinations, which are 
unimportant for the qualitative behavior of the spectrum (see corresponding figures in 



Shibanov et al. 



1992 



Suleimanov et al. 



20091 ). Besides, the bright regions at magnetic poles 



of a star in our models considered below have magnetic fields close to normal. 

The emergent spectra of two thin magnetized model atmospheres (i.e. a hydrogen layer 
above a condensed surface) with parameters typical for the XDINSs are shown in Fig. HJ 
We display the Eddington flux 



He = - / Ie{cos a) cos a sin a da, 
2 Jo 



(13) 



where Je is a specific intensity. Later we also use the astrophysical flux = 4:He- The 
often used physical flux is defined as J-'e = vtFe = in He- For blackbody radiation it is 
-^E = -^E = Be- In the models shown in Fig. 2 the blackbody radiation is used as the inner 
boundary condition for the radiation transfer equation. There are absorption features in 
the spectra, which arise due to the proton cyclotron absorption and bound-bound and 
bound-free transitions in hydrogen atoms. It is necessary to remark that the proton 
cyclotron feature is weakened by vacuum polarization and partial mode conversion, but 
nevertheless remains pr ominent in the spectra, as can be seen in our figures below (see also 



Suleimanov et al. 



20091 ). Equivalent widths (EWs) of these features are ~ 25 eV and ~ 90 
eV for the models with surface densities S = 1 and 10 g cm~^, respectively. The absorption 
feature in the spectrum of a semi-infinite atmosphere is even less prominent in corapariso n 



Suleimanov et al. 



20091) 



to the absorption feature in the model with E = 10 g cm~^ (see 
and its absorption feature has a smaller EW (~ 50 eV). The EW of the absorption feature 
slightly depends on the magnetic field strength and the relation between proton cyclotron 
energy and kT^fi- For example, the EW of the absorption feature in the spectrum of the 
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model with B = 10^^ G and S = 10 g cm ^ and the same T^s is close to 80 eV (see Fig. H]). 

The spectra of these models can be fitted approximately by diluted blackbody spectra 
Be with a color temperature Tc slightly larger than the effective temperature Teg and two 
Gaussian absorptions: 

Fe = D BeUc Teff) exp(-ri) exp(-r2), (14) 

where fc = Tc/T^s is the hardness factor, D is the dilution factor, which can differ from the 
usual Later we will present the dilution factor in the form D = f^~^ to distinguish it 
from the usual value. The optical depths of the absorptions at a given photon energy E are 



(E — Ei^2 



n, = r° exp -^^-4^ . (15) 



2<2 

Here ^^e optical depths at the centers of the lines at energy £'1^2 and o"i,2 are the line 

widths. The values of these fitting parameters for these two spectra are given in the caption 
of Fig. |2l The first line corresponds to the proton cyclotron line in the model atmosphere 
spe ctra, and the second one co r responds to the boun d-bound transitions in hydrogen atoms 



(cf. IPotekhin fc Ghabriej 120041 : IPotekhin et al.ll2004J ). 



It is clear that the blackbody radiation model is not an accurate inner boundary 
condition for the thin model atmospheres above the condensed surfaces, and the real 
radiation properties of the surface have to be used for precise modeling. Below we use our 
approximations (|9]) - f|TT]) for the inner boundary conditions of radiation transfer equation 
for both radiation modes, extraordinary (X) and ordinary (O): 

I^ISV.E) = ^5£;(T(w))(l - 2R) + 2R I^lJ-ri.E), (16) 

/,Kr/,E) = ii?^(T(w)). (17) 

where rj = cos a, and / are the specific intensities of the extraordinary and 
ordinary modes. In future work we plan to replace this approximation by accurately 
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Fig. 2. — Emergent spectra of two thin model atmospheres above a blackbody condensed 
surface (sohd curves) with Tgg = 1.2 x 10^ K and B = 4x 10^^ G together with fitting spectra 
(dashed curves). The fitting spectra are blackbody spectra with the color temperature 
Tc = fcTcs and two Gaussian lines. Top panel: The spectrum of the model with surface 
density S = 10 g cm^^. Parameters of the fitting spectrum: fc = 1.2, D = 1.14^^, center of 
the first line (proton cyclotron) Ei = 0.25 keV, optical depth ti = 2.6, and width (Xi = 6 eV, 
and for the second line (H bound-bound transitions) E2 = 0.3 keV, T2 = 1.1, (72 = 40 eV. 
Bottom panel: The spectrum of the model with surface density S = 1 g cm~^. Parameters of 
the fitting spectrum: /c = 1.12, D = 1.12~^, line parameters are the same as above, except 

T2 = 0.08. 
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computed radiation properties of the condensed surfaces. But we expect that the present 
approximation is quahtatively correct. 

We have computed thin hydrogen model atmospheres above condensed iron surfaces 
using these new inner boundary conditions. The results are presented in Figs. [3]-[8l Model 
atmospheres are optically thick in the 0-mode in the observed part of the X-ray band (0.1 
- 5 keV), and almost everywhere optically thin in the X-mode. Therefore, the total spectra 
are defined mostly by the X-mode spectra, because the 0-mode spectrum arises in the 
cooler layers of the atmosphere. 

The emergent spectra of the thin magnetized model atmospheres above the condensed 
surface are harder (/c ~ 1.5 - 1.7) in comparison to the models above a blackbody, because 
the temperatures at the inner boundaries of these models are higher (see bottom panels in 
Figs.|3]-|l]). 

Due to the condensed surface emission properties, there appear broad absorption 
features between £'c,i(Fe) and approximately Eq, especially prominent at Ec^\. The 
dependence of optical depths of thin atmospheres in both modes on the photon energy 
are shown in Fig. [51 The atmospheres are optically thin in X-mode at almost all energies, 
except the region around the proton cyclotron energy. The total spectra are dominated by 
the X-mode (see Fig. [6]), therefore, they are close to the condensed iron spectra in X-mode, 
especially in the Wien tail. But this is not the case for the band between i?c,i(Fe) and 
Eq. In this region the 0-mode flux can contribute a significant part to the total flux (see 
Fig. E]). At energies in the vicinity of -E'c,i(H) the opacity in the X-mode increases, therefore, 
the optical depth of atmosphere and the emergent flux increase, too. As a result a complex 
absorption feature arises with a total EW ~ 400 eV for the model with B = 10^^ G and S 
= 10 g cm-2, EW ^ 370 eV for 5 = 4 x lO^^ G and S = 1 g cm'^, and EW ^ 260 eV 
for B = A X 10^^ G and S = 10 g cm~^. It is interesting to remark that the EWs of these 
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0.1 1 
E (keV) 




Column density (g cm"^) 

Fig. 3. — Emergent spectra (top panel) and temperature structures (bottom panel) of thin 
partially ionized hydrogen atmospheres {T^s = 1.2 x 10^ K, B = A x 10^^ G) above a 
condensed iron surface with S = 10 g cm"^ (solid curves) and S = 1 g cm~^ (dashed 
curves). An inner boundary condition for the radiation transfer equation, corresponding 
to the condensed surface radiation properties (see text), was used. The blackbody spectra 
with T = Tgfj (dotted curve) and the diluted blackbody spectra (/c = 1.5, D = 1.3^'^), 
which describe the high-energy part of the spectrum of the atmosphere with S = 1 g cm~^ 
(dash-dotted curve) are shown for comparison in the top panel. In the bottom panel the 
temperature structure of the model with S = 10 g cm~^ using blackbody radiation as the 
inner boundary condition is also shown (dotted curve). 
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Fig. 4. — Effect of the inner boundary condition of the radiation transfer equation on emer- 
gent spectrum (top panel) and temperature structure (bottom panel) of a thin partially 
ionized hydrogen atmosphere (Teg = 1.2 x 10^ K, B = 10^^ G) with S = 10 g cm~^, above 
the condensed iron surface (solid curves) and blackbody surface (dashed curves). For com- 
parison the blackbody spectrum with T = Teg (dotted curve) and the diluted blackbody 
spectrum (/c = 1.6, D = 1.4~^), with a Gaussian line representing the proton cyclotron line 
(-El = 0.63 keV, ri = 8.5, ui = 11 eV), and a half-Gaussian line {E2 = 0.3 keV, ri = 2.5, (Xi 
= 150 eV) describing the absorption feature at energies higher than (dash-dotted curve) 
are shown. That last curve fits the spectrum of the atmosphere above the condensed iron 
surface. 
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Fig. 5. — Optical thickness of a thin partially ionized hydrogen atmospheres with Tgs = 
1.2 X 10^ K above the condensed iron surface in X-mode (solid curves) and 0-mode (dashed 
curves). Top panel: B = 4 x 10^^ G and S = 1 and 10 g cm~^. Bottom panel: B = 10^^ G 
and S = 10 g cm~^. 
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Fig. 6. — Total emergent spectra (solid curves) together with spectra in the X-mode (dashed 
curves) and in the 0-mode (dotted curves) of a thin partially ionized hydrogen atmosphere 
with = 1.2 X 10® K, S = 10 g cm~^ above the condensed iron surface with B = 4 x 10^^ 
G (top panel), and B = 10^^ G (bottom panel). 
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E (keV) 

Fig. 7. — Emergent specific intensity spectra of two model atmosplieres above a condensed 
iron surface witli Tcs = 1.2 x 10^ K, S = 10 g cm"^, and B = A x 10^^ G (top panel) and 
B = 10^^ G (bottom panel). Spectra at different angles between radiation propagation and 
surface normal are shown: 1° (solid curves), 30° (dashed curves), 60° (dotted curves) and 
75° (dash-dotted curves). 
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Fig. 8. — Angular distributions of the emergent radiation at photon energy 1 keV. Top panel: 
Models with T^s = 1.2 x 10^ B = Ax 10^^ G above the condensed iron surface and S = 
10 g cm~^ (solid curve), S = 1 g cm~^ (dashed curve), S = 10 g cm~^ above a blackbody 
(dotted curve), and a semi- infinite atmosphere (dash-dotted curve). Bottom panel: Models 
with Tgfj = 1.2 X 10^, B = 10^^ G and above the condensed iron surface with S = 10 g 
cm~^ (solid curve), S = 1 g cm~^ (dashed curve), S = 10 g cm~^ above a blackbody (dotted 
curve), and a semi-infinite atmosphere (dash-dotted curve). 
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features are larger than the sums of the EWs of the cyclotron line in the atmosphere spectra 
and the EWs of the absorption feature in the condensed surface spectra (see Eq. ([7])). The 
reason is the following. Between E^^i and Eq a naked condensed surface radiates mainly in 
the 0-mode, which is approximately half of the flux in comparison to other energies. In 
contrast, a thin atmosphere above it is optically thick in the 0-mode, and the emergent 
flux forms in the upper atmosphere layers, which have significantly lower temperatures than 
the condensed surface. Therefore, the emergent flux in this band is significantly smaller 
compared to the naked condensed surface. This is particularly noticeable in the model with 
B = 10^^ G. 

Apart from the spectral energy distribution of the emergent flux, the angular 
distribution is also important. Examples for various models of magnetized atmospheres 
are shown in Figs. [7] - [HI There is a narrow spike at the surface normal, but the total 
energy radiated in this spike is relatively small and we don't consider these spikes below. 
Angul ar distributions of semi- infinite magnetized atmospheres have maxima at a ~ 40° 



60° (iPavlov et al. 



I994J ). An example of such a distribution is presented in the bottom 



panel of Fig. El It can be fitted by the function 

/£;(«) = Ee X 0(a), 

where 



[a 



3.6 cosa — 0.5 cos^ a — 2.9 cos^ a. (19) 



The angular distributions of the specific intensities of the thin magnetized atmospheres 
above a blackbody are close to isotropic (Fig. [8]). The corresponding distributions for the 
atmospheres above condensed surfaces are also almost isotropic ones (see FigJH]), but have 
some interesting feature at energies between Eq and 4£'c. The location of the high energy 
boundary of the broad depression depends on a (see Eq. f ll2p ). Therefore, at a given photon 
energy there is an angle a corresponding to this location. At this angle the specific intensity 
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angular distribution has a jump. At smaller angles the distribution corresponds to the 
usual angular distribution for the model atmosphere, but at the larger angles the specific 
intensities are significantly smaller (see Fig. |H]). This jump is more prominent for the model 
with B = 4 X 10^^ G and it is also well visible in the spectra of the specific intensities 
for this model (the top panel of Fig. [7]). In the angular distributions of the intensities of 
the model with B = 10^'' G this jump is less significant (Figs. [7]and|H]), because for this 
model Eq ~ Ec^\{H) and the X-mode opacity is more important, together with the mode 
conversion. The features of the magnetized model atmospheres are very unusual and they 
have to be studied in more detail with a more accurate treatment of the radiation properties 
of the condensed surfaces. 

From the presented examples of the partially ionized hydrogen magnetized model 
atmospheres we can conclude, that they display complex absorption features with EW 
< 100 eV in the spectra. The naked condensed surfaces can have broad absorption features 
with slightly larger EWs (> 120 eV), but the absorption structures in the spectra of thin 
atmospheres above condensed surfaces can have significantly larger EWs, up to 400 eV. 



3. Calculation method for the integral neutron star spectrum 

We consider a slowly rotating spherical isolated neutron star with given compactness 
M/R as a model for XDINSs spectra and light curves fitting, taking into account 
distributions of the local temperatures and magnetic field strengths over the stellar surface. 
Local spectra are blackbody-like spectra with corresponding temperature and absorption 
lines, whose central energies depend on the local magnetic field. Gravitational redshift and 
the light bending are taken into account. 

The geometry of our model is shown in Fig. |9l The magnetic axis is inclined to the 
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rotation axis by the angle ^b, and i is the inchnation angle of the rotation axis relative to 
the line of sight. The phase angle is the angle between the plane containing the rotation 
axis and the line of sight and the plane containing the magnetic and rotation axes. One 
magnetic pole can be shifted relative to the antipode of the other pole by a small angle k. 
We assume that the brightest areas of the NS surface are located at the magnetic poles. 
Therefore, temperature and magnetic field distributions near the poles are most important. 

XDINSs pulse fractions are significant, and a classical temperature distribution, which 



arises due to a global (co re) dipole r aagnetic field o 



observed pulse fractions (jPage 



1995 



Geppert et al 



a N S, is not sufficien t to re produce 



2006f ). 



Schwope et al 



fl2005h showed 



that the RBS 1223 light curves can be fitted by two models, namely two bright isothermal 
spots with different temperatures , or two narrow peak e d temperature distrib utions (with 



different parameters) taken from 



Geppert et al 



too 



1- 



Geppert et al 



(I2OO4I ) considered 



the thermal transport in the NS crust for a poloidal dipole magnetic field, concentrated in 
the crust only. In this case the magnetic field lines have large inclinations to the surface 
normal in a wider band around the magnetic equator than in the global dipole field case. 
The thermal conductivity across magnetic field lines is smaller than along them, therefore 
the bright spots around magnetic poles are brighter and smaller than they would be for 
the dipole distribution. Another possibility is a strong toroidal m a gnetic field in the NS 



Perez- Azorm et al 



(l2006al ) using various 



crust. Such kind of models were computed by 
approximations. They obtained an analytical approximation of the surface temperature 
distribution for the crust magnetic field with a strong toroidal component calculated for the 
force-free conditions. We use this approximation for the temperature distribution in our 
model. 



In a well known approximation by iGreenstein fc Hartkd (119831 ) the temperature 
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Fig. 9. — Geometry of the model. The opposite bright spot is indicated by the dashed circle. 
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distribution is: 



COS $ H sm $ 



(20) 



where i^x and i^y are thermal conductivities across and along the magnetic field, 
respectively, and Tp is the temperature at the magnetic pole. Magnetic fields of the isolated 
NSs can be strong enough (~ 10^'^ G) so that the ratio K\ / K \\ is small. In this case the 



temperature distribution can be approximated as (jPerez-Azorm et al. 



2006bh : 



T^ = Tp^cos^$ + T^i,. (21) 
In the force-free approximation a simple relation between t he magnetic colatitude 9 an d the 



magnetic field inclination in the NS crust was suggested by 



Perez- Azorm et al. 



2006ah : 



2 4 cos^ 6 

"""^ " (l + /x2i?2) + (3_^2^2)cos2^' ( ) 

where is a typical toroidal magnetic field length scale, and /ii? is an approximate ratio of 
the toroidal component of the magnetic field to the poloidal one. The case /x=0 corresponds 
to a core dipole magnetic field. Although this magnetic field model was derived for the 
crust, we use Eq. fl22l) also in the atmosphere, because it allows us to investigate the possible 
effects of magnetic field and temperature inhomogeneity by changing a single parameter /i. 
Finally, we use the following temperature distributions around magnetic poles: 

rr.A _ ^4 ^"^^^ I 7^4 

cos^ U + ai^2 sm 6 

where ai,2 = (1 + A*i2-^^)/4 and Tpi^2 are the temperature distribution parameters. This 
can describe various temperature distributions (see Fig. [10]), from strongly peaked (a 1) 
to the classical dipolar (a = 1/4) and homogeneous (a = 0). 

Another possibility is to consider two isothermal bright spots around the magnetic 
poles with fixed temperatures and angular sizes. 
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Fig. 10. — Temperature distributions at various parameters a. Only the first term on the 
right-hand side of Eq. ( 123|) is used for these calculations. 
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The surface magnetic field distribution, which corresponds to the magnetic field 
inclination in the force-free approximation model (Eq. [2^ is described as follows: 

B = 5pi,2-\/cos2^ + ai,2 sin^^, (24) 

where -Bpi,2 are the magnetic field strengths at the poles. 

Temperature and magnetic field distributions are calculated separately around each 
magnetic pole, therefore, the distributions are not smooth at the magnetic equator area if 
one of the poles is shifted. But this small mismatch at the equator is unimportant, because 
this area is dim. 

The local spectrum of the NS surface in the case of magnetized atmospheres at 
given temperature and magnetic field is taken as a diluted blackbody spectrum with two 
additional absorption features (see also Eqs. (fT4l) . f|T5|) and (fT8|) ) 

Je(«) = D BMcT) 0(«) exp(-ri) exp(-r2), (25) 

where the optical depth Xj of the absorption line i at given photon energy E is defined 
by Eq. f|T5|) . and r°2 5 ^1,2 are considered as model parameters identical at all NS surface 
elements. The photon energy of the line center Ei depends on the local magnetic field 
strength B. In particular, for the proton cyclotron line we have 

El = 0.0635 Bn keV. (26) 

For the description of the local spectra of the model atmospheres above a condensed iron 
surface, two additional absorption features can be considered. The first feature is centered 
at the iron ion cyclotron energy according to Eq. ([T]). This line is asymmetric, i.e., r2 > 
only for E > E2, in agreement with the calculations (see Figs. OandS]). The second line 
is centered between Ei and Eq, and its position depends also on the local magnetic field 
strength. 
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We use three angular distributions (f){a) of the specific intensity /E(a). The simplest 
one is the isotropic function (f){a)=l, which corresponds to the thin model atmosphere 
above a blackbody. For the semi-infinite model atmosphere the angular distribution f lT^ 
are used. For the thin atmosphere models above a condensed iron surface the isotropic 
angular distribution is used, except for the energy band between Eq and 4i?c- Here the 
angular distribution has a jump at some cos a (see Fig. |H]), and this case has to be studied 
separately. 

For the description of the naked condensed surface spectra the relations ([H]) - f[T^ are 
used. 

The total spectrum of the NS is cal culated in the spherical coord inate system connected 



with the rotation axis (see, for example 

?2 />27r 



Poutanen fc Gierlihski 



2003 



and references therein) 



/e = / dip (1 + 2;)"^ Je' (a) COS a sin 7^7, (27) 
Jo Jo 

where d is the NS distance and 7 is the angle between the radius vector at a given point and 
the rotation axis. Here the observed and the emitted photon energies are connected by the 
relation E = E'{1 + z)~^, and z = {1 — Rs/R)"^^"^ — 1 is the gravitational r edshift. Ligh : 



bending in the gravitational field is accounted by the approximate relation (IBeloborodov 



20021 ) 



Rs cos ip 
"R ^ {1 + z 

where Rs = 2GM/c^ is the Schwarzschild radius, and 



COS a ~ — + y:;— ; ^, {Zi 



cosV' = cosicos7 + sinisin7cos(y9. (29) 

The integration in Eq. ( 127|) is restricted to the visible NS surface area with cos a > 0. 

Phase resolved spectra can be calculated for various phase angles, and the corresponding 
light curves can also be calculated for various energy bands. The mean spectrum is summed 
from separate phase spectra and then divided by the number of phases. 
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Using this code we calculated two test light curves from slowly rotating neutron stars 



(see F ig. ITT]) . The first one corresponds to the light curve presented by 



Poutanen &: Gierlihski 



( 120031 ) in their Fig. 4a (dashed curve), a nd the second one c orresponds to Class III light 



curves in the classification performed by 



Beloborodovl (|2002|) (see his Fig. 4). 



4. Results of modeling 

In this section we present some examples of model light curves and integral emergent 
spectra for magnetized neutron stars with various surface temperature distributions. We 
choose symmetric models with equal temperatures and magnetic field strengths at both 
magnetic poles, and equivalent temperature and magnetic field distributions in both 
hemispheres. We consider a neutron star model with polar temperature kTp = 0.15 keV, 
polar magnetic field strength i?p = 6 x 10^'^ G and a compactness corresponding to the 
gravitational redshift z = 0.2. Three temperature distributions are used. Two of them are 
described by Eq. ( 123|) with a = 0.25 and a = 60. In both distributions Tmin equals 0.316Tp. 
The results are almost the same for any smaller value of Tmin. The third temperature 
distribution is a uniform surface temperature kT = 0.01 keV plus two uniform bright spots 
around the magnetic poles with the same temperatures Tp and angular radii 9sp = 5°. Most 
of the models are calculated for angles 1 = 6-^ = 90°, which gives the largest change of 
the visible local parameters during a rotation period. For simplicity we take dilution D 
and the hardness fc factors equal to 1. Most of our results will not be changed by this, in 
particular, it will not change the normalized spectra and EWs of absorption lines, if fc and 
D are constant over the neutron star surface. This approximation will change the absolute 
value of the flux only. Therefore, it is important for the ratio of X-ray to optical fluxes, 
because fc and D can be different for X-ray and optical bands (see Sect. 4.3). Therefore, in 
this section we will consider the temperature distribution over the neutron star as a color 
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Fig. 11. — Test li g ht cu rves calculated with the same p arameters as used by 
Poutanen fc Gierlihskil (120031 ) (top panel) and iBeloborodovl (120021 ) (bottom panel) . Fi is 
the maximum possible flux that is observed when one of the spots is viewed face-on. 



temperature distribution instead of the effective temperature distribution. 



4.1. Averaged spectra 

The main goal of our modehng is the study of possible absorption features in the 
magnetized neutron star thermal spectra. To this aim the spectra of the models with the 
three temperature distributions, as described above, and the three local spectral energy- 
distributions, as described in Section 2, are computed. In all the cases we consider the color 
temperature distribution over neutron star surface. 



4- 1.1. Blackbody plus Gaussian line 



The first spectral model is a simple isotropic blackbody spectrum with a temperature 
equal to the local temperature of a given model and one Gaussian spectral line. The line is 
centered at the proton cyclotron energy (see Eq. (126!) ) and has parameters ri = 3, (Xi = 30 
eV, hence EW ~ 100 eV. This model roughly represents the shape of the spectrum of a thin 
magnetized partially ionized hydrogen atmosphere above the blackbody condensed surface 
(see Fig. [2]). The same model can represent the shape of the spectrum of the semi- infinite 
atmosphere at the maximum of the flux, but, of course, with larger fc. Model spectra 
of semi-infinite magnetic atmospheres have hard tails in comparison with a blackbody. 
But to first approximation these spectra at the maximum flux can be described as 



diluted blackbody spectra (see, for example. 



Suleimanov &: Werner 



20071 for non-magnetic 



atmospheres). Moreover, the model spectra with vacuum polarization and partial mode 
conversion have smalle r deviations from black body than the models, where these effects are 



ignored (see Fig. 10 in 



Suleimanov et al. 



20091 ). We mainly investigate here the absorption 



features, and ignore the harder spectral tail in the real model atmosphere spectra. 
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Fig. 12. — Phase averaged photon spectra for the neutron star model with kTpi = kTp2 = 
0.15 keV, Bpi = = 6 X lO^^ G, z = 0.2, i = Ob = 90° and different temperature 
distributions across the surface. Top panel: the temperature distribution Eq. ( |23l) with 
ai = a2 = 0.25 (dotted curve), with ai = a2 = 60 (dashed curve), and two uniform bright 
spots with T = Tp and angular size Ogp = 5° (solid curve). The local spectra are isotropic 
blackbody spectra with an absorption line (ri = 3, cti = 30 eV) centered according to 
Eq. (126 p . The spectrum of the model with a = 0.25 and a weaker absorption line (ri = 2.6, 
cxi = 6 eV) is also shown (dash-dotted curve). The thin solid curves are the continuum 
spectrum of a model with a = 60 and a single blackbody spectrum with a Gaussian line 
{El = 0.25 keV, ri = 0.85, ai = 52 eV) for the strong-line a = 0.25 model. Bottom panel: 
The photon spectra for two temperature distributions (a = 0.25 and two bright spots) under 
various inclination angles to the line of sight: i = 0° (solid curves), i = 45° (dashed curves), 
and i = 90° (dotted curves). 
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Integral photon spectra for all the temperature distributions are shown in the top panel 
of Fig. [121 The overall shapes of the spectra can be described by blackbody spectra with 
kT ^ 0.114 keV for the model with a = 0.25, kT ^ 0.086 keV for the model with a = 60, 
and kT 0.125 keV for the model with two bright spots. Therefore, these spectra can be 
described in first approximation by single blackbody spectra with a Gaussian line. The 
example of such an approximation is shown in Fig. [T^l The model spectrum with a = 0.25 
is well represented by the blackbody spectrum with kT ^ 0.114 keV and a Gaussian line 
with El = 0.25 keV, ti = 0.85, cti = 52 eV. The integral spectra have a soft excess and an 
asymmetry of the line in comparison with single blackbody spectra. These features increase 
if a increases (see also Fig. fTTjl . 

The spectrum of the model with a = 60 is the softest one due to a relatively large 
contribution of the cooler parts of the neutron star. The spectrum of the two spots model 
has a relatively narrow absorption line, because most radiation comes from the pole regions 
where the magnetic field strength is nearly constant. The spectrum of the model with 
a = 0.25 shows a wide asymmetric line, shifted to lower energy. The large width of the line 
is a consequence of the smooth temperature distribution, and, accordingly, a large variety 
of magnetic field strengths, which effectively contribute to the spectrum. The absorption 
feature in the spectrum of the model with a = 60 is very wide and difficult to see. It 
can be recognized by a comparison with the spectrum without the line (the thin curve 
in Fig. [121 top panel). In this model the magnetic field changes significantly, from Bp at 
the poles up to 7. 7 Bp at the equator and the absorption line is strongly smoothed. An 
observable absorption feature appears only in models with a < 10 - 50 (depending on the 
EW of the local absorption line). Therefore, strong toroidal fields on the surface of XDINSs 
(corresponding to a > 50) are incompatible with the observed absorption lines. 

The equivalent widths of these absorption features range from 65 eV (a = 60 model) 
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to 85 eV (the other models). Therefore, the magnetic field distribution does not result in 
an increase of the EW of the absorption line in the integral spectrum. For illustration the 
integral spectrum of the model with a = 0.25 and a line with ti = 2.6, ai = 6 eV and EW 
~ 19 eV is also shown in the upper panel of Fig. [121 The line has the significant width but 
a small depth with EW ^17 eV. 

The integral spectra of the model with a = 0.25 and the two-spots model at different 
inclination angles of the magnetic axis and rotation axis (which coincide), to the line of sight 
{i = 0°, 45° and 90°) are shown in the bottom panel of Fig. [T21 These spectra represent 
phase resolved spectra at various phases of the stellar rotation for the case 6*3 = 90°. The 
spectra of the two-spots model only show the fiux level variation, without a change of the 
line shape. In contrast, the spectra of the models with a = 0.25 show line shape variations. 
At larger i the line center is shifted to smaller energy and becomes wider, because a more 
significant part of the observed radiation comes from the magnetic equator regions. The 
dependence of the absorption line EW on the inclination angle is insignificant. 

4.1.2. Naked condensed iron surface 

The second model used for the local spectrum is the spectrum of the naked condensed 
iron surface, provided by Eqs. ([8]) - f|T2|) . The computed spectra for the neutron star 
models with the same parameters and temperature distributions, as in the case of the first 
local model, are shown in the top panel of Fig. [131 These integral spectra are close to 
blackbody spectra with similar color temperatures, as in the first case (see Sect. 4.1.1). 
The absorption features are very wide and shallow and have EWs equal to ~ 110 eV for the 
two-spots model, ~ 140 eV for the model with a = 0.25, and ~ 95 eV for the model with 
a = 60. In the last case the line is again not observable, and the EW is obtained with a 
lower accuracy because of uncertainties in the continuum definition. 
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Fig. 13. — Top panel: Pliase averaged plioton spectra of tlie same neutron star models as 
in Fig. [121 but liere tlie local spectra are from a naked iron surface. The continuum spectra 
without absorption features are shown by the thin solid curves for the two-spots model and 
the model with a = 0.25. Bottom panel: Phase averaged photon spectra of the same neutron 
star models as in Fig. [121 but here the local spectra are isotropic blackbodies with absorption 
features like in the spectra of the thin hydrogen model atmosphere above the condensed iron 
surface (see text). The continuum spectrum without absorption feature is shown by the thin 
solid curve for the model with a = 60. 
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It is possible to compare approximate calculations of neutron star spectra with naked 
co ndensed iron surface with_the accurate calculations, performed by the method described 

( 2005! ) ■ An example of this comparison is shown in Fig.[TU The 



m 



van Adelsberg et al. 



qualitative agreement is satisfactory. The accurate spectrum is (as expected) more smooth 
at Eq, and the flux at energies less than Ec^\ is larger. The EWs of the absorption features 
are very similar, 190 eV for the accurate spectrum and ~ 180 eV for the approximate 
spectrum. 



4- 1.3. Thin atmosphere above condensed iron surface 

The shape of the spectrum of a thin hydrogen atmosphere above a condensed iron 
surface can be roughly fitted by a blackbody spectrum with two absorption features with 
Gaussian profiles, which describe the proton cyclotron line and one half of the wide line 
at i?c,i(Fe). The half line means that the line optical depth is equal zero at energies 
E < i?c,i(Fe), and calculated in the usual way (Eq. (fT5|) ) a.t E > i?c,i(Fe). For our illustrative 
calculations we take ri = 8.5, ai = 11 eV, and T2 = 2.5, (J2 = 150 eV with a total EW ~ 
230 eV (see also Fig.Hj). The Gaussian lines with these Tj and well fit the absorption 
features in the spectrum of the thin hydrogen atmosphere above the condensed iron surface 
with Teff = 1.2 X 10^ K, 5 = 10^^ G, and S = 10 g cm'^ (Fig. i]). The computed spectra 
for the used temperature distributions are shown in the bottom panel of Fig. [131 Again, 
the temperature of the neutron star surface means the color temperature. The hard tails 
of the integral spectra can be fitted by the blackbody spectra with approximately the same 
temperatures as those obtained for the first local spectral model. The absorption feature 
has a prominent two-component structure for the two-spots model, and a smoothed shape 
for the model with a = 0.25. In the spectrum of the model with a = 60 the absorption line 
is again very smooth. The EWs of the features in spectra of the two-spot model and the 
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Fig. 14. — Phase averaged photon spectra of the neutron star models with parameters 
kTpi = kTp2 = 0.15 keV, Bpi = Bp2 = 6 x 10^^ G, i = 9b = 0° and a = 0.25. The local 
spectra are from a n aked iron surface, computed accurately using the code presented by 



van Adelsberg et al. 



( l2005l ) (thick solid surve), and using approximate Eqs. ([H]) - (IT^ (thick 
dashed curve). The adopted continuum spectra without absorption features are shown by 
the corresponding thin curves. The EWs of the absorption features are rather similar, ^190 



and 180 eV. 
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model with a = 0.25 are ~ 210 eV, and the EW of absorption feature in the spectrum of 
the model with a = 60 is ~ 150 eV. 



4.2. Light curves 



Light curves calculated for illustrative purposes in a wide energy band (0.02-2 keV) are 
shown in Fig. [151 The light curves are computed for the orthogonal rotator {i = 6b = 90°), 
the three temperature distributions, and two angular distributions of the emergent flux: 
the isotropic one and the distribution corresponding to a semi-infinite magnetized hydrogen 
model atmosphere as described by Eq. (fTOl) . The light curves, computed with the isotropi c 
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angu lar distribution, agree with the well known fact (see, e.g. 
20061 ) that the temperature distribution corresponding to a core dipole magnetic field 
(a = 0.25) cannot provide a pulsed fraction PF = (F^ax — -^min)/(-^max + -^min) larger 
than several percent. In order to ac count for larger PF values, more peaked temperature 



distributions have to be considered (ISchwope et al. 



20051 ). Computed light curves have PF 



~ 5% for the temperature distribution with a = 0.25, ~ 10% for a = 60, and ~ 25% for the 
two bright spots. 

Using the angular distribution corresponding to the magnetized model atmosphere, we 
obtain completely different pulse profiles with four maxima instead of two for the isotropic 
angular distributions (bottom panel in Fig. [T^ . The model with a = 0.25 shows additional 
maxima at the phases of the minima in the light curve computed with isotropic angular 
distribution. The models with the narrow peaked temperature distributions (two spots 
and a = 60) show deep minima at the phases of the maxima in the light curves computed 
with isotropic angular distribution. Each minimum in these light curves becomes replaced 
by two maxima with a shallow minimum between them. The pulsed fractions of the light 
curves for the models with a = 0.25 and a = 60 are the same as in the isotropic distribution 
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Fig. 15. — Examples of light curves for the neutron star model with fcTpi = kTp2 = 0.15 
keV, Bpi = Bp2 = 6 X 10^'^ G, z = 0.2, i = 6-q = 90° and different temperature distributions 
across the surface: the temperature distribution fl23l) with ai = a2 = 0.25 (dotted curves), 
with ai = a2 = 60 (dashed curves), and two uniform bright spots with T = Tp and angular 
size ^sp = 5° (solid curves). Two different angular distributions of the emergent radiation 
are used: the isotropic distribution (top panel) and the semi-infinite magnetic atmosphere 
distribution Eq. ( IT9l) (bottom panel). 
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case, but the pulsed fraction of the hght curve for the model with two bright spots increases 
up to 60%. A qualitatively similar Class III light curve for a slowly ro tating neutron star 



with magnetic atmosphere model radiation was obtained by 



Hd (120071 ) (see his Fig. 2) 



We also investigate the influence of the absorption features on light curves and pulsed 
fractions. This influence on the light curves in a wide energy band is very small. But it is 
significant on the light curves and pulsed fractions in the narrow energy bands at a vicinity 
of the absorption feature energies. For illustration we calculate the dependences of PF on 
the photon energy for the models with a smooth temperature distribution (a=0.25) and 
the different local spectra (see Fig. [TS]). The corresponding averaged spectra are presented 
in the same figure for illustrative purpose. The common characteristic of the all presented 
dependences is the decreasing of the pulsed fraction at the absorption feature energies. 

In the top panel of Fig. [16] the local spectrum is approximated by the blackbody with 
and without an absorption line, and the changing in the pulsed fraction at the absorption 
line energies very well seen. The PF for the bolometric light curve (~ 5%) is determined 
by the soft part of the spectrum, where most of the photons come. The PF can be more 
significant for the higher energy bands. We also calculate the averaged spectrum and the 
PF dependence on energy using the same model but a more peaked angular distribution of 
the local emergent radiation for an electron scattering atmosphere: 

0(a) = 0.4215 + 0.86775 cos a. (30) 

The averaged spectrum is changed insignificantly, but the pulsed fraction is increased by a 
factor of two (the top panel in Fig. [T6|) . We present this result in order to demonstrate that 
a pencil-like peaked angular distribution of the emergent radiation can give a similar effect 
like a peaked temperature distribution. 

The dependence of the pulsed fraction on energy for the models with local spectra 
corresponding to a naked iron surface and a thin atmosphere above the condensed iron 
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surface is similar (the bottom panel in Fig. [T6|) . The PF at energies E > Eq is larger due 
to peaked angular distributions (see Eqs. ([H]) - ffT2]) and Fig. [H]). In our calculations we 
approximate the angular distribution of the thin atmosphere above the condensed surface 
emergent radiation by a simple step function in the Eq - AEq energy range: 



\ ^ ~ a^cosa^ 

(pia) = it cos a > cosOc, (ol) 

1 — cos 



and 

0(a) = a^, if cosa < cosoc, (32) 

where 

COS = 77 1 , (33) 



and is a parameter, which was taken 0.2 here. At other energies (j){a)=l is used. 

The pulsed fraction at these energies can be even more significant if slightly different 
temperatures at the poles are used. Here we do not take into account an azimuthal 
dependence of the angular distribution of the magnetized semi-infinite atmospheres with 
inclined magnetic field. Therefore, the light curves can be even more complicated. But the 
observed shapes of most XDINS light curves are close to a simple sine, therefore we expect 
that the angular distribution of the emergent radiation from XDINS is close to isotropic. 



4.3. Fluxes in the optical band 

Here we compare the fluxes of our computed models in the optical band with the 
fluxes of the blackbody extrapolations of the X-ray model spectra. It is well known that 
the observed optical fluxes of XDINSs exceed the blackbody extrapolations many times 
(see Introduction), therefore, this comparison can give additional information on the 
temperature distribution and the nature of local spectra. The comparisons are presented 
in Fig.[T7l The models with a smooth temperature distribution (a=0.25) show the small 
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Fig. 16. — Energy dependence of pulsed fractions (symbols and left vertical axes) together 
with averaged photon spectra (lines and right vertical axes) for the neutron star models with 
the same parameters as in Fig. [12] and a=0.25. Top panel: The models with pure blackbody 
local spectra (open circles and dotted curve), blackbody local spectra with an absorption line 
(r = 3.0, (7=30 eV) centered according to Eq. fl26l) . and the isotropic angular distribution 
(filled circles and solid curve), and blackbody local spectra with the absorption line and the 
electron scattering atmosphere angular distribution (triangles and dashed curve). Bottom 
panel: The model of a naked condensed iron surface (open circles and dotted curve), a thin 
atmosphere above condensed iron surface spectra with corresponding angular distribution 
(filled circles and solid curve), and the same local spectra, but with slightly different pole 
temperatures: Tp i = 0.15 keV, Tp 2 = 0.14 keV (triangles and dashed curve). The position 
of Eq, which corresponds to the magnetic pole, is shown by vertical dotted hne. 
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excesses (factor ^ 1.3 - 1.5) for all local model spectra. The narrow peaked temperature 
distributions show more significant excesses (up to a factor of ~ 10), which depend on the 
values of a and the adopted temperature for the cool part of the neutron star surface. 

There is an additional uncertainty connected with dilution D and hardness /c factors. 
In our magnetized neutron star model spectra we represented the local model atmosphere 
spectra by blackbody spectra with D = /c = 1 for simplicity. But in the real model 
atmosphere spectra these factors can dif fer from 1 and, moreover, they can be diffe rent for 



the X-ray and the optical band spectra (iPavlov et al. 



1996 



Suleimanov et al. 



20091). 



Conclusions 



With this paper we have started our study of the XDINS surface radiation properties. 
We have considered a few theoretical models of a highly magnetized {B > 10^^ G) neutron 
star surface: partially ionized hydrogen atmospheres, semi-infinite and thin; condensed iron 
surfaces; and thin partially ionized hydrogen atmospheres above a condensed iron surface. 



by 



To describe the e merg ent spectra of the model atmospheres, we have used the results 



Suleimanov et al. 



(120091 ) and demonstrated that the equivalent widths of the absorption 



features do not exceed 100 eV. We have presented a simple analy tical approximatiori of the 



eme rgent spectra of the c o ndens ed surfaces, previously studied by 



Perez- Azorm et al. 



and 



van Adelsberg et al. 



(120051) 



and evaluated the equivalent widths of the absorption 



features in these spectra (120 - 190 eV for XDINS magnetic fields, see Eq. ([7])). 

Using this approximation, we have studied thin partially ionized hydrogen atmospheres 
above a condensed iron surface. The radiation properties of the condensed surface are 
used for the inner boundary condition for the model atmospheres. Some examples of such 
models are presented. We have demonstrated that the interaction of the condensed surface 
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E (keV) 

Fig. 17. — Comparison of the integral averaged spectra of the magnetized neutron star 
models (solid and dashed curves) with blackbody spectra, which fit the model spectra at 
E > 0.5 keV (dotted curves). The differences at the optical band are also shown. Top panel: 
The models with blackbody local spectra and the Gauss line (r = 3.0, cr=30 eV), centered 
according to Eq. (126 p . for three temperature distributions: a=0.25 (upper curves), a=60 
(middle curves), and two-spot model (lower curves). The parameters of the distributions are 
the same as in Fig. [151 The spectrum of two-spot model with the temperature A;T=0.03 keV 
of the cold surface is also shown (dashed curve). Bottom panel: The models with smooth 
temperature distribution (a=0.25) and various local model spectra: naked iron spectrum 
(solid curves) and the thin model atmosphere above condensed iron surface spectra (dashed 
curve). Parameters of the local spectra are the same as in Fig. [131 
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radiation and the radiation transfer in the atmosphere leads to strongly increased equivalent 
widths of the absorption features, up to 300 - 450 eV. Another interesting property found 
is the prominently peaked angular distribution of the radiation at energies between Eq and 
4^0- 

All of these different local models used to compute integral emergent spectra of 
neutron stars with parameters close to XDINS ones. We presented the modeling code for 
light curves and integral spectra. We have considered a magnetic field distribution with 
a possible toroidal component and two different models for the temperature distribution 
on the neutron star surface. The first model is a simple uniform relatively cold surface 
with two uniform hot spots with given angular s izes at the magnetic poles. The second 



model was taken from 



Perez- Azorm et al. 



(l2006al ) and describes a continuous temperature 



distribution, peaked at the magnetic poles. An additional parameter a (connected with a 
relative contribution of the toroidal component of the magnetic field) controls the width of 
the temperature profile peaks. 

Light curve modeling confirms that a classical core dipole temperature distribution 
(a = 0.25 in our model) with local blackbody spectra is not sufficient to explain the observed 
pulsed fractions in XDINSs. Therefore, more peaked temperature distributions or more 
peaked angular distributions of the specific intensity are needed. We have demonstrated 
that the isotropic or peaked angular distribution of the emergent radiation, as well as 
angular distributions in the models of thin atmospheres, better explains the observed 
pulse profiles in comparison with the magnetized semi-infinite model atmosphere angular 
distribution. 

We also investigated dependence of the pulsed fractions on energy. The pulsed fraction 
increases with photon energy due to a strong dependence of the thermal flux on the 
temperature in the Wien part of the spectrum. The pulsed fraction decreases at absorption 
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feature energies. The emergent flux angular distribution of a thin atmosphere above the 
condensed iron surface and a naked condensed surface is peaked to the surface normal at 
an energy range between Eq and 4:Ec, therefore, the pulsed fraction is increased there. It is 
also increased if the pole temperatures are different. 

We computed some examples of integral emergent spectra of hot magnetized neutron 
star models for various temperature distributions and all considered surface models. We 
show that the absorption features in the integral spectra become wider, but do not show an 
increase of the equivalent widths due to the variation of the magnetic field strength over the 
neutron star surface. Therefore, the observed EWs in XDINS spectra have to be similar to 
the local EWs on the neutron star surfaces. Most of the XDINSs have observed absorption 
lines with EWs < 100 eV, which can be explained by magnetized model atmospheres or 
naked condensed surfaces, but the EW of the absorption line in the spectrum of RBS 1223 
is about 200 eV, and an atmosphere above the condensed iron surface provides a better 
explanation. 

The models with a strong toroidal magnetic field on the neutron star surface (a > 
10 - 50) display a too wide smoothed absorption feature, which cannot be detected by 
observations. Therefore, a strong toroidal magnetic field component on XDINSs surfaces 
can be excluded. At the same time the large pulsed fraction in some XDINSs hint to the 
existence of small hot regions on the neutron star surfaces. Th erefore, the model w ith a 



poloidal magnetic field concentrated in the neutron star crust (iGeppert et al. 



2004h and 



models w ith a strong toroida. 



example. 



Aguilera et al. 



component in the crust, which vanishes on the surface (for 
20081 ). can be applicable. On the other hand, similar result can 
be reached if the peaked angular distribution of the emergent spectra is considered. For 
example, the thin hydrogen atmosphere above the condensed iron surface with a smooth 
temperature distribution over the neutron star surface (a=0.25) and slightly different pole 
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temperatures can provide the pulsed fraction observed in RBS 1223. 

In the present work the radiation properties of the magnetized neutron star condensed 
surfaces are considered in a rough approximation only. More accurate investigations of the 
magnetized model atmospheres above the condensed surfaces are necessary in the future. 
However, our present model can be used for fitting of observed spectra of XDINSs to obtain 
basic information about radiation properties of magnetized neutron star surfaces. 
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